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ABSTRACT: Complex systems such as glasses, gels, granular materials, and 
systems far from equilibrium exhibit violation of the ergodic hypothesis (EH) 
and of the fluctuation-dissipation theorem (FDT). Recent investigations in 
systems with memory [1] have established a hierarchical connection between 
mixing, the EH and the FDT. They have shown that a failure of the mixing 
condition (MC) will lead to the subsequent failures of the EH and of the FDT. 
Another important point is that such violations are not limited to complex 
systems: simple systems may also display this feature. Results from such sys- 
tems are analytical and obviously easier to understand than those obtained 
in complex structures, where a large number of competing phenomena are 
present. In this work, we review some important requirements for the validity 
of the FDT and its connection with mixing, the EH and anomalous diffusion 
in one-dimensional systems. We show that when the FDT fails, an out-of- 
equilibrium system relaxes to an effective temperature different from that of 
the heat reservoir. This effective temperature is a signature of metastability 
found in many complex systems such as spin-glasses and granular materials. 

1 Introduction 

Since its formulation by Boltzmann [2], the EH has called the attention of 
mathematical physicists and chemists. In the last century, a branch of the 
mathematics dedicated to its study has been developed. However, most of its 
results are accessible only to the specialist. On the other hand, the FDT has 
played a central role [3,4] in nonequilibrium statistical mechanics in the linear 
response regime (NESML). It gained such importance that Kubo proposed a 
complete formulation of the NESML based on it [4] . Since the FDT is directly 
related to relaxation processes, its more empirical character has caught the 
attention of experimentalists and most of the discussion about its validity 
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has remained in the hands of theoretical physicists and chemists, instead of 
mathematicians. 

A necessary requirement for the vahdity of the FDT is that the time- 
dependent dynamical variables be well defined at equilibrium. The presence 
of nonlinear effects or far from equilibrium dynamics may lead to situations 
where the FDT does not hold [1,5], the aging process in spin-glass systems 
being a good example [6-11]. 

Most of experimental situations in which the EH and the FDT arc violated 
happen in complex structures. Nevertheless, we show here simple situations 
where those violations appear. Their explicit condition and simplicity allows 
the judgment of non-specialists in the subject. This is important because it 
opens the possibility that more complex structures can be investigated on a 
solid basis. 

This work is organized as follows: In this section, we shall define the EH, 
the MC, and other main concepts and ideas to be discussed throughout this 
work. In Section 2, we outline some historical achievements in the study of 
diffusion and introduce the FDT. Wc discuss reaction rates in Sec. 3, be- 
fore introducing, in Sec. 4, the concept of memory, i.e., we discuss a system 
governed by a Generalized Langevin's Equation (GLE) and we show how to 
obtain anomalous diffusion. Next, wc discuss random walks, fractional deriva- 
tives and their connection to the GLE in Sec. 5. We then continue by defining 
in a clear manner the noise in Sec. 6, and its connection with memory and cor- 
relation functions. After that, wc discuss reversibility in Sec. 7. In Section 8, 
we discuss the main issue of this review, which is the interconnection between 
the MC, the EH, and the FDT, and we show under what conditions they 
fail. Examples of such violations are given in Sec. 9, where we study ballistic 
motion. In Section 10, we introduce some speculative topics on the forefront 
of physical research; the "skeptic reader" can skip it. Finally, we introduce a 
conjecture in Sec. 11 and conclude the paper in Sec. 12. 

Let us start by considering the evolution of a dynamical stochastic variable 
A{t) see Fig. 1. The variable could be either at equilibrium, Fig. 1(a), or 
approaching it, Fig. 1(6). The ensemble average {G{A)) of any function G{A) 
is defined as 



where = fc^T is the inverse temperature, E{A) are the energies, and the 
integration is performed over all the accessible states of the phase space i?. 
From that, it is possible to define a correlation function as 



For an exponential decay of the correlation function, it is possible to associate 
a relaxation time r, which is larger than the typical time for a fluctuation At. 
One can also define a time average as 




(1) 



CA{t) = {A{t)AiO)). 



(2) 
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Fig. 1. Evolution of a dynamical stochastic variable. 
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For Z\t <C r T, the average will produce the continuous line in Fig. 1, i.e., 
it will wash out the fluctuations which we measure with a sensitive probe. For 
times T ^ T, the EH reads 



A{t) = {A{t)). (4) 

In simple words: given enough time, the system will reach every accessible 
state, and a time average will be equal to an ensemble average. The proof 
of the EH has been the Holy Grail of statistical mechanics. One expects it 
will hold for macroscopic systems at equilibrium, Fig. 1(a), and for small 
deviations from equilibrium; on the other hand, it will not hold for curve 
(6), although it is expected that the system in situation (6) will be driven 
to equilibrium for longer times. However, a general proof of the EH is still 
missing. The concept of "far from equilibrium" is itself sometimes misleading, 
since it depends not only on the initial conditions, but also on the possible 
trajectories the system may follow [1, 12]. The way a system approaches an 
equilibrium is crucial for these definitions. 

The "mixing property of a physical system" or mixing condition (MC) can 
be stated as 

lim R{t) = 0, (5) 

t — *oo 

where we have defined the normahzed quantity R{t) = CA{t)/CA{0). The MC 
tells us that after a long time t^T,we do not expect that A{t) will remember 
its initial value A{0). 
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2 Diffusion 

Diffusion is one of the most fundamental mechanisms for transport of energy, 
mass and information; it is a main process for a system to reach uniformity and 
equihbrium, and has therefore been the focus of extensive research in many 
different disciplines of natural science. For almost two hundred years, it has 
caught the attention of the scientific community. The famous observations by 
Robert Brown [13, 14] of the erratic trajectories of pollen opened a new world 
for experimental and theoretical studies in what was named Brownian motion. 
As a biologist. Brown first assumed he had discovered the basic essence of 
life, an idea to be expected from a man dedicated to biology. Nevertheless, 
he reproduced the experiments with non-organic material and observed the 
same erratic motion. Brown then concluded that this was due to the motion of 
matter. Considering his aim, the second conclusion is not only more difficult; 
it is a highly advanced and honest conclusion. Unfortunately, most books are 
unfair with Brown, in not mentioning his subsequent experiments. 

At the centennial celebration of the Einstein miraculous year, one could 
be easily driven to the conclusion that most of the diffusive phenomena are 
well understood today. However, if we ask simple questions such as "How do 
spin waves diffuse in a Heisenberg system with correlated disorder?", "How 
do electrons behave in an irregular lattice?", or "How does a ratchet device 
work?" , it takes a short time to realize that these unanswered problems are 
related to diffusion. 

For instance, when we flip a spin in the ground state of a ferromagnetic 
chain, the principle of equal a priori probability for the accessible states tells 
us that somehow the energy due to this disturbance will not remain localized 
in a single state. However, it does not say whether or not the system will 
support a spin wave, whether the wave propagates, and how it propagates in 
the affirmative case. To answer this kind of question, we usually need to go on 
into specific calculations. Our aim is to try to understand the general character 
of diffusion and, hopefully, to classify it prior to extensive calculations. 

At the beginning of the twentieth century, much important research was 
dedicated to the irregular motion of microscopic particles dispersed in a fluid, 
namely Brownian motion. Despite being irregular, the motion reveals some 
regularity when analyzed statistically. The main observed quantity was the 
mean square displacement {x^{t)) of the particles, which evolves linearly with 
time. 

Einstein's basic idea was to explain the Brownian motion by going be- 
yond thermodynamics and into kinetic theory. He considered single spherical 
particles of mass m and radius a suspended in a liquid of viscosity r], and 
obtained [15] 

lim {x'^{t)) = 2Dt. (6) 

t — ^oo 

As usual, by infinite time, wc mean a time larger than the maximum relaxation 
time of the process. The diffusion constant D was given by 
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D = = = 

where R is the gas constant, Na the Avogadro number, fj, the mobiUty, and 7 
the friction the particle feels in the fluid. 

Equation (6) was a major achievement; the linear relation with time was 
confirmed and an expression for the diffusion constant was obtained. The 
last form for D in Eq. (7) establishes a connection to the mobility, which is 
fundamental for the study of conductivity and transport. It was then possible 
to check the theory with the data available at the time. If one knows D, it 
is possible to make an estimation of the Avogadro number Na- Future works 
helped establish the Boltzmann constant, ks = R/Na, as a. new fundamental 
constant. Moreover, estimation of the size of sugar molecules dissolved in water 
became possible. Besides that, the frequency-dependent diffusion constant 
D{u)) can be directly associated with the conductivity by the relation [16] 

= ^Dico), (8) 

where e is the carrier charge and n the carrier density. 

The series of Einstein articles about diffusion [15] together with the work of 
Smoluchowski [17] paved the way for the modern theory of Brownian motion. 

The next step forward was taken by Langevin. In order to describe the 
motion of a particle immersed in a fluid, in 1908, Langevin [18] proposed the 
equation 

= -m-/v{t) + f{t). (9) 

where f{t) is a stochastic force subject to the conditions {f{t)v{0)) = 0, 
{f{t)) = and {f{t)f{t')) = A6{t - t'). If we solve Eq. (9) and by using 
the equipartition theorem we impose — > 00)) = fc^T/m, we obtain the 
proportionality constant A = 2mkBTj and write 

(/(*)/(*')> = 2mkBT^5{t - t'). (10) 

This last relation is known as the fluctuation-dissipation theorem (FDT). 
Although the Einstein diffusion constant contains implicitly the relation be- 
tween fluctuation and dissipation, in Langevin's formulation, it acquires the 
importance of a basic theorem. 

From Eq. (9) and the above conditions, one obtains the velocity-velocity 
correlation function 

Cy{t) = {v{t + t')v{t')) - (fesT/m) exp(-7i). (11) 

The correlation function Cv{t), or R{t) = exp(— 7i), will satisfy Eq. (5), the 
MC, with a relaxation time r = 7~^. This exponential decay from the initial 
conditions is the expected form for the MC. Now, it is possible to obtain the 
mean square displacement as 



6 Mendeli H. Vainstein, Ismael V. L. Costa, and F. A. Oliveira 



{x''{t > r)) = /* dt' f dt"{v{t')v{t")) = 2Dt, 



(12) 



Jo Jo 

where, in the last step, we used the Kubo formula 




oo 



Cv{t)dt. 



(13) 



The Kubo formula, together with Eq. (11), reproduces Einstein's results. 

There are many reasons why one should always look back to Langcvin's 
work, the first one being that it focuses attention on the motion of a particle, 
which is very intuitive for any physicist. Second, it combines the old Newto- 
nian deterministic approach with the new "uncertainty" of the stochastic force 
f{t). The breaking of the atomic forces in two parts: a fast changing force f{t) 
with time scale At, and the slow friction force with time scale r, introduces 
a large simplification which facilitates understanding and computer simula- 
tion. Consequently, the use of the Langevin equation, and of its generalization 
(Sec. 4), is still very active, having been applied successfully to the study of 
many different systems such as the dynamics of dipolar [19] and polymeric 
chains [20-23], metallic liquids [24], Lcnnard Jones liquids [25], diffusion in 
periodic potentials [26], ratchet devices [27,28], and synchronization [29], only 
to name a few. Finally, it established explicitly for the first time the connec- 
tion between fluctuation and dissipation, the FDT, which remains a major 
theorem of statistical mechanics. 

The Langevin equation, however, presents some limitations: (a) It is a 
classical formalism; (b) It has uncorrelated noise with only two time scales 
At and r, whereas a complex system has in general many time scales; (c) We 
cannot make any predictions for times shorter than At; (d) It predicts only 
normal diffusion. 

A quantum formidation of the FDT has been put forward by Callen and 
Welton [30]. Following their work, much research has been done in the field, 
with many attempts at generalization [6,31]. We shall focus our attention on 
the Kubo FDT [4,32,33], or the so-called second FDT, since it is more useful 
to the study of diffusion; see Sec. 4. 

3 Reaction Rates 

Diffusion may be considered the simplest problem of nonequilibrium statistical 
mechanics; however, if one considers a particle moving in an irregular media, 
it is quite probable that the particle will be affected by some potential and, 
in moving from point A to point C, will have to cross a potential barrier at 
point b, see Fig. 2. This is one of the oldest problem in statistical mechanics 
and was first mentioned by Arrhcnius [34] . Arrhcnius proposed that the rate 
of particles crossing the potential barrier will be given by 



k = i'oexp{-/3Eb) 



(14) 
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Fig. 2. Potential with two minima. 

where uq is the attempt frequency and Eb is the barrier height. Arrhenius as- 
sociated the attempt frequency to the vibration frequency at the well bottom. 

The issue stood still until Kramers [35] re-addressed it with a more accu- 
rate analysis, in such way that the problem is known today as the Kramers 
problem. He gave an expression for the attempt frequency, which depends 
on the curvature at the top of the potential barrier and on the friction. The 
Kramers theory of reaction rates has many applications in biology, chemistry 
and physics. For a review, sec [36]. 

In a recent work, Sancho et al. [26] studied the diffusion of particles in a 
periodic potential. They simulated the diffusion using Langevin's equation and 
used Kramers reaction rate to analyze the jumping between successive wells. 
A good agreement between theory and the simulation was obtained. This is a 
nice exhibition of the relation between transport processes and reaction rates. 

For systems with memory [37 39], the Kramers problem suffers the same 
restrictions as those of diffusion; worse, some of the restrictions are not even as 
explicit as in diffusion [40]. In short, the reaction rate theory will be enriched 
by a better understanding of diffusion. 

4 Complex Systems have Memory 

If one asks "What makes a system complex?", a few concepts will come to 
mind: a large number of degrees of freedom, nonlinearity and memory. How- 
ever, if we ask for a definition of memory, only a few answers will be precise. 
We shall use here the concept of memory introduced by Mori in his seminal 
paper [37], where he used a method of projection operators which has many 
advantages: it allows the treatment of quantum systems, it is not empirical, it 
has time correlation, and it is a non-Markovian formulation with an explicit 
definition of memory. 
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Forty years after Mori's work, many ftindamental concepts and methods 
have been developed [3, 4, 37, 41 47], which allowed a generalization of the 
Langevin formalism and the elimination of most of its limitations (see items 
(a) to (d) at the end of Sec. 2). The new formalism gives origin to a Generalized 
Langevin Equation (GLE) of the form 



dA{t) 
dt 



= - f n{t-t')A{t')dt' + F{t), 

Jo 



(15) 



where F{t) is a stochastic noise subject to the conditions {F{t)) = 0, 
{F{t)A{0)) = and 

CF{t) = (F(i)F(O)) = {A%,n{t). (16) 

Equation (16) is the Kubo FDT [4,32,33], being a generalization of Eq. (10). 
The memory, n{t), arises here explicitly and, in principle, it allows us to study 
a large number of correlated processes. Notice that now many time scales are 
possible within 7T(t); this is a natural condition for complexity. An equation 
for CA{t), or for the renormalized correlation function R{t) is given by 



dR{t) 
dt 



f n{t-t')R{t')dt', (17) 
Jq 



where we have used the conditions {A{0)F{t)) = 0. The Laplace transform of 
this equation yields 

z + n{z) 

Prom here on, we shall use the tilde to indicate Laplace transforms. 
Let us now define the variable y{t) as 

y{t)= f A{t')dt', (19) 
Jo 

with asymptotic behavior 

lim (2/2 (t)) ~ r. (20) 

For normal diffusion, a = 1; we have subdiffusion for a < 1, and superdiffusion 
for a > 1. Notice that if A{t) is the momentum of a particle, then y{t)/m is 
its position. 

It is very simple to show that this new formalism allows both normal and 
anomalous diffusion. Consider two examples: first, take 

n{t) = 2-f5{t). (21) 

With this short range memory, we return to the normal Langevin's equation, 
Eq. (9), and obtain a = 1, i.e., normal diffusion. Second, consider an extremely 
long memory 
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n{t) = iol (22) 

which gives a force of the form —mtOQ-y, i.e., an harmonic oscillator which does 
not exhibit diffusion at all, a = 0. Those are artificial but simple examples of 
how the memory determines diffusion. 

Recently, Morgado et al. [48] obtained a general classification for anoma- 
lous diffusion. They used the time-dependent diffusion function 

D{t)= I Cy{t')dt', (23) 



^0 

lim Dit) = lim zD(z) = lim n~^(z) = lim ^ . (24) 

t— >oo 2— >0 z— ►0 t— »oo 77(1/^) 

Here, we have used the final value theorem [49], and Eq. (18). Note that 
the long range time behavior, i.e. the dynamics, is dominated by the small 
frequencies. We expect that z ^ faster than r{z). Consider now that 

n{z ^ 0) ~ z". (25) 

With t ~ z~^, we get 

a = u+l. (26) 

This result is fundamental for the classification of diffusion. Let us take 

the previous examples: first, n{t) = 2jS{t) with Laplace transform 7T(z) = 7; 
here v = 0, a =_1, and the diffusion is normal as expected. Second, from 
Eq. (22), we get n{z) = oJq/z, v = -1, and a = 0. 

We are now in condition to discuss mixing. Consider a system governed 
by a GLE. The asymptotic behavior 

lim n{t) = lim zn{z), (27) 

lim R(t) = lim zR(z) = lim ^ , (28) 

t^oo z^o z^o 2 + n{z) 

and condition (25) show that mixing exists only for — 1 < < 1. In other 
words, the mixing condition is fulfilled between the limits of ballistic motion 
and harmonic oscillator motion. Close to the limits, one can expect problems. 
Notice that for the memory given in Eq. (22), we have the exact solution 

R{t) = cos{ujot), (29) 

which obviously does not fulfil the MC, Eq. (5). Notice also that if Eq. (13) 
converges, then it is always possible to have an associated "friction constant" 
7, even for correlated systems 



7 



/•OO 

-1= / R{t)dt. (30) 
Jo 
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This definition can be plugged into the diffusion formula (7). This explains 

why normal diffusion is so widely found in nature, even for processes we 
clearly know are strongly correlated [48]. For a long time, the only evidence 
of anomalous diffusion had been for subdiffusion processes due to trapping 
mechanisms [16] and hierarchical lattices [50]. However, supcrdiffusivc motion 
has recently been studied both theoretically [51,52] and experimentally [53-58] 
(see Sec. 9). We also expect that research in the new nanoclay technology will 
produce anomalous diffusion [59]. Recent works on chain dynamics [19,22,60] 
show that the system dynamics may build up a memory. This "casual" result 
is explored in the conjecture described in Sec. 11. 

5 Random Walk 

The study of random walks is somehow older than statistical mechanics and it 
has produced many alternative ways to describe diffusive processes. Even be- 
fore Einstein, the mathematical works of Bernoulli opened up the possibility 
of understanding fluctuations; for example one can arrive at Eq. (6) by consid- 
ering the famous drunk man problem, a starting point in many undergraduate 
texts [61]. Besides that, the famous law of large numbers, or N~^^'^, for the 
relative standard deviation of a variable suggested that if atoms exist they 
must be very small, with N being very big, otherwise thermodynamics would 
not make sense. However, the famous expression for diffusion Eq. (7) was 
not obtained before Einstein; it has the basic information one needs to know 
about the nature of the process. That is a main difference between mathemat- 
ics and physics. In this sense, the Einstein diffusion constant was a stunning 
achievement. We shall call the attention to the famous Chandrasekhar review 
on stochastic process [62], which was an up-to-date article until the origin of 
the Mori formalism. It still remains as a clear and concise review. 

By the end of the nineteen century, the works of Lord Rayleigh [63] on 
random flights allowed to understand scattering in random directions, one of 
the most famous application of which was to explain the blue color of the sky, 
due to light scattering by impurities in the atmosphere. A modern formulation 
shows the latter approach is incorrect and that light scattering is associated 
with fluctuations in the dielectric constant [64]. After the creation of the laser, 
light scattering became itself a large field of research. [65,66] 

Let us consider here a simple random walk analogy to Langevin's work. 
Consider a set of A''o particles with initial velocities equal to zero. The parti- 
cles are subject to a random force ±/o at each time interval At with equal 
probability. After N time steps, the average velocity will be {v{t)) = and 
the average squared velocity 

{v'{t)) = NflAe/m\ (31) 

Wc can see from the previous expression that the kinetic energy grows linearly 
with time, t = N At, i.e., the random force acts as a pump of energy. Such 
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a simple model does not represent a physical process, since it contradicts 
the kinetic theory, — ^ oo)) ~ kBT/m. To make it more realistic, we 
add a dissipative force —m'yv{t), and impose the balance of energy to obtain 
/o = 2mkBTj/At. Notice that if we define the Dirac delta function as the 
limit At ^ oi S{t) = I /At for -At/2 < t < At/2, and 0, otherwise, we 
recover the FDT, Eq. (10). The FDT is nothing more than a detailed balance 
condition; it is a guarantee that the dispersed particles will reach thermal 
equilibrium after a reasonable time. 

In the same way, there are alternative ways to describe anomalous diffu- 
sion besides the GLE. One proposal is to use the continuous random walk, 
which can be mathematically described by fractional derivatives [67-69] . The 
fractional derivative of a function f{x,t) can be defined as 

1 d f'^^, fix,t') 



oDj-°'f(x,t) = —-— dt' /^ ' / , (32) 
* ■'^ ' ' r{a)dtJo {t-t'y-"' ^ ' 

where r{x) is the Gamma function. Equation (32) is a natural generalization 
of the derivative of a complex variable using the residue theorem. The nonlo- 
cal character of the fractional derivative is the same as that of the memory. 
Therefore, it is quite natural that they yield the same results as the GLE. 
Indeed, it is possible to obtain a fractional Fokker-Planck equation (FFPE) 
of the form 



d V'{x) d 



f{x,t), (33) 



dx mrja dx"^ 

to address the problem of subdiffusion [68] and to obtain a relation similar to 
Eq. (26). 

The generalization for superdiffusion was recently discussed in a few arti- 
cles ( [69], and references therein). Fractional derivatives are a very compact 
way to obtain results; however, there are some points one should bear in mind. 

First, we assume a priori a fractional geometry and as a result we obtain 
fractional diffusion; that appears to be a circular argument. From the GLE, 
that comes naturally from the memory, or from the noise; see Eq. (36). 

Second, if x{t) gives the time evolution for a particle position, we know 
precisely what oD'//x{t) means only for /it = 1, 2, . . . However, no one has an 
idea of what it means for /i = 0.51, or for any non-integer value. 

Finally, natural solutions for the fractional derivatives are the Levy func- 
tions. Unfortunately, these yield an infinite mean square displacement, which 
is not a good physical result. The Levy distributions, (^^(.t), are very popular 
because they fulfill the generalized central limit theorem (GCLT) [70] 

cl>^{x) = j (t>^{x - x')(t>^{x')dx' . (34) 

Moreover, recently [71,72], it has been shown that the GCLT represents the 

first uncorrelated term in a renormalization process. Correlations, such as 
those one expects to find in anomalous diffusion, will cause the deviation 
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of the studied variables from the GCLT. Recently, Figueiredo et al. [73-76], 
have proposed some theorems on the limit sum of stochastic variables without 
making the classical assumptions of the GCLT. They developed a general 
formalism to explain the non- convergence (or the slow convergence) to the 
Gaussian distribution. With this, they have explained the origin of the self- 
similar property that appears in real economics time series data. They have 
also explained how autocorrelations (linear and nonlinear) can be considered 
as a soiirce of truncated Levy flights. The asymmetry they have found in their 
distributions are similar to those found in relaxation in supercooled liquids 
and in the height distributions in the etching of a crystalline solid [77-79]. The 
experiments of Monte et al. [80,81] show both asymmetric and superdiffusive 
behavior. 

6 Noise 

A fundamental aspect of stochastic processes is the noise. A stochastic gener- 
alized force F{t) can be decomposed into a set of harmonic oscillators of the 
form 



Here, < 0(a;) < 2-k is the random phase. From the FDT, it follows that 



where p(w) is the noise density of states (NDS). Now, we take the Laplace 
transform of Eq. (36) to obtain 



This is another relevant result. The friction is equal to the noise density of 
states at the origin. This shows how the lower modes determine the type 
of diffusion. A system which has a finite friction presents normal diffusion, 
since its NDS is finite at the origin. Subdiffusion will have an infinite friction 
and an infinite NDS. Finally, superdiffusion has a null friction. To obtain 
super diffusion for a null friction is a very intuitive and appealing concept. A 
null NDS tells us that the lower modes do not relax and the process has "less 
interference" , or, in Langevin's language, "weak collisions" . 
Consider now the colored noise 



F{t) = {A^)l'g^ / p{iof'^cos{ujt + (t>{'^))<Lj 



(35) 





(37) 




(38) 



Here, lod is a Dcbye cutoff frequency. This kind of noise has been used by 
Caldeira and Leggett in quantum dissipative systems [82]. If we plug this 
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noise into Eq. (36), take its Laplace transform, and then the Hmit of small z, 
we obtain the exponent [83] from Eq. (25) 

.= 1^ '^f^<l (39) 
II , otherwise. 

For most of the cases, the exponent of the NDS for low frequencies will 
be the same as that of the Laplace transform of the memory for small z. 
Equation (39) shows that a < 2 and, consequently, the motion is limited by 
the ballistic motion. Ballistic motion appears to be a limit of this kind of 
GLE, see Sec. 8. Notice as well that for z/ = 0, we get 7 = 70, from Eq. (37). 

We shall consider here another possibility. Let the noise be 

I U , otherwise. 

For coi = 0, we have the Debye density of states for a thermal noise com- 
posed of acoustic phonons. Thus, for wi = we have normal diffusion and for 
any lui =^ we have superdiffusion. This NDS is the difference between two 
Ornstein-Uhlenbeck processes and is a simple way to produce ballistic diffu- 
sion [1,48]. Since there is a window, < w < wi, where there is no fluctuation 
of the modes, this introduces a very practical mechanism to control simu- 
lations. This kind of noise seems more appropriate to describe real ballistic 
propagation [84] than Eq. (38), see Sec. 9. 



7 Reversibility and Correlation functions 

As early as 1876, Loschmidt called attention to the reversibility paradox [85]. 
His paradox states that all molecular processes must be reversible, since there 
is a symmetry between past and future t ^ —t in the laws of physics. Con- 
sequently, statistical mechanics must be reversible, in apparent contradiction 
with thermodynamics, where certain processes are irreversible. At that time, 
reversible physics was composed of classical mechanics and electrodynamics. 
This paradox, together with the dynamical problems of Liouville, Zermello, 
and Poincare are central in the work of Boltzmann. Those lead to the Boltz- 
mann equation, to the H theorem [86], and to the studies of the Poincare 
recurrence theorems. 

Again, our aim here is not to go into extensive mathematical proof; rather 
wc! focus on the correlation function. In the definition of the correlation func- 
tion we have used 

(41) 

and 

R{-t) = R{t). (42) 
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The first relation relies on our basic knowledge of the temporal invariance 
of physical laws. However, if the evaluation is made far from equilibrium, 
the correlation function R may depend both on ti and t2, and not on their 
difference. 

The second relation is time reversal, which can be easily understood for 
a classical variable where A{t) and ^(0) commute. Given a string of values 
A{ti) i = 1,2,..., Nint, for large Nmt, one can obtain the relations given by 
Eqs. (41) and (42). For quantum systems, the reader is recommended the 
review of Balucani et al. [87] 

A great achievement in the discussion of time reversal symmetry for macro- 
scopic systems was made by Onsagcr. He considered a solid subject to a gen- 
eral field E. In the linear regime, the field induces a generalized current density 
J of the form 

Ji = o-i,kEk- (43) 

The susceptibility a satisfies 

0'i,k = IJ'i,k<^k,i- (44) 

Here, ^i^i = 1 for the diagonal terms. The off-diagonal terms are ^i^k = 1 for 
variables which do not change sign under time reversal, such as the electric 
field, and = — 1 for variables which do change, such as the magnetic field. 
Equation (44) is the Onsager reciprocal relation. 

A natural generalization of susceptibility is the correlation function or the 
response function, sometimes called a Green function. Consider, for example, 
that one applies a perturbation P{xi) at the position xi and wants to know 
the disturbance S{x2) at X2- In the linear regime, we get 

Si{x2) = j Gi^k{x2,xi)Pk{xi)dxi. (45) 

For systems with translational invariance, we expect that 

G{xi,X2)=G{xi-X2). (46) 

When the translational invariance is broken, due to the existence of surfaces 
such as in a film [88,89], defects, or topological disorder [90], the response 
becomes a function of both variables xi and X2- However, Oliveira [88] has 
given a proof based on time reversal symmetry that 

Gi,k {xi,X2) = Hk,iGk,i {x2 , xi ) (47) 

The tensor /Xjj here is the same as the one in Eq. (44), i.e., we lose the space 
invariance, but, on the other hand, wc gain a useful space exchange symmetry. 
This result is general and has many applications. For example, it was used to 
explain certain asymmetries found in light scattering [88]. 

The usual Onsager reciprocal relations are actually limited to systems 
asymptotically close to equilibrium. For example, they apply at the level of 
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the Navier-Stokes equations for a simple fluid, but fail for the higher-order cor- 
rections to those equations, as pointed out by McLcannan [91,92]. Dufty and 
Rubi [92] generalized Mcleannan's work to nonequilibrium stationary states. 

Many correlation functions of the form {A{xi,ti)A{x2-t2)) have properties 
similar to Eq. (46). Notice that even in nonlinear systems, sometimes it is 
possible to make some general statements. For example, in the growth process, 
the height of a surface h{x, t) is a function of the position x and of the time t. 
The main studied quantity is the roughness [93], defined by the mean square 
fluctuation 

Ah'^{x, t) = {{h{x, t) - {h{x, t))f). (48) 
The roughness satisfles the scaling laws 

bAh^{bx, ¥-t) = Ah^ix, t). (49) 

where 6 is a number and z is the growth exponent. Notice that this relation 
holds only statistically. Many symmetries or scaling in the correlation function 
hold in situations where nothing can be said for a single process. 

We now return to the correlation functions of the GLE. Note that the 
memory is an even function of t, independent of the NDS (see Eq. (36)). The 
analytical extension of the Laplace transform of an even function is an odd 
function, n[—z) = —IJ{z). Consequently from Eq. (18), R{—z) = —R{z), and, 
by a converse argument, R{t) is an even function [83]. This is in agreement 
with the results by Lee for Hamiltonian systems [94]. Notice also that Eq. (17) 
requires the derivative of R{t) to be null at the origin. Since both memory 
and R{t) are even, they can be written as 

oo 

n{t) = Y,bnt''', (50) 

n=0 

and 

oo 

i?(i) = ^a„t2". (51) 

n=0 

Unfortunately, a large number of works has been presented in the litera- 
ture where the correlation function does not satisfy these requirements. We 
shall not comment further on these works here: some may be useful approx- 
imations, others represent artificial solutions. The reader should be cautious 
in identifying them. 

Exponentials, stretched exponentials, and power laws are examples of 
asymptotic behavior that can be obtained from more complex even func- 
tions [83], but obviously they do not fulfil Eq. (17). 

Determination of the coeSicients in Eqs. (50) and (51) can be done for 
every specific noise. For short times, those equations yield R{t) = cos{LUot), 
where loq = 11(0). For broadband noise, the asymptotic times yield expo- 
nential decay for normal diffusion. For anomalous diffusion, the behavior is of 
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a stretched exponential followed by inverse power-law. For short band noise a 
very rich oscillatory behavior may be found [83] . 

A large and growing literature in which non-exponential behavior has been 
observed for correlation functions can be found in the following articles: in 
glasses and supercooled liquids [95], frustrated lattice gases [96], liquid crys- 
tals [97-99], plasmas [53], hydrated proteins [100], growth [101] and disordered 
vortex lattice in superconductors [102]. 

8 Mixing, Ergodicity, and the Fluctuation-Dissipation 
Theorem 

In this section, we arrive at the central point of this work by showing that 
the EH, Eq. (4), the MC, Eq. (5), and the second FDT, Eq. (16) are strongly 
connected in the GLE. Consequently, one could expect that the violation of 
one of these conditions could lead to the violation of the others. However, we 
will show that there is a hierarchy among the three concepts, in such a way 
that some may be violated, while others are not. Most of the systems that 
present violation of the FDT are complex, such as supercooled organic liquids, 
some algebraic maps, evolutionary models, and the eternal spin-glass problem. 
We try to show here a minimal condition for violation of that hierarchy. 

We may expect from Eqs. (15) and (16) that a system will be driven to 
an equilibrium state, i.e. 

lim A^) = (52) 

which can be identified with the EH. We shall see, however, that this is not 
always the case for superdiffusive dynamics. 

Note that the Laplace transform of Eq. (15) suggests a solution of the form 

A{t) = A{0)R{t) + I R{t- t')F{t')dt', (53) 
Jo 

where we have an ensemble of initial A{0). Squaring Eq. (53) and taking the 
ensemble average, we obtain for the asymptotic behavior [1] 

{A^{t ^ 00)) = {A^),, + R\t ^ oo)[(^2(o)^ _ ^^2^^^]^ (54) 

This simple result leads to very important consequences. First, the system will 
reach full equilibrium only if the MC, Eq. (5), holds. Second, the EH holds if 
the MC holds. Finally the FDT will hold only if the EH holds. Consequently, 
the FDT is the end validation of the sequence: Mixing ^ Ergodic Hypothesis 
Fluctuation-Dissipation Theorem. Observe that if the MC is violated, then 
the final value of Eq. (54) will depend on the initial conditions. That is just 
the essence of the MC. 
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At this point, we shall call again the attention to Lee's work in ergodic- 
ity [45]. Unlike any other previous attempt at establishing the validity of the 
Boltzmann EH, his work approaches the time average directly and explicitly, 
which was made possible by his recurrence relation method [43]. Moreover, it 
has been demonstrated in several exact solvable models when the hypothesis 
is valid. When it is not valid, it is shown the reason why the hypothesis fails. 

If the mean square value of A can be associated with a given temper- 
ature by the equipartition theorem, we have (A^(0)) ^ To for the initial 
temperature, T for the reservoir temperature, and Te// for the final effective 
temperature. Equation (54) becomes 

T,ff=T + R^t^(x^)[To-T]. (55) 

From Eq. (28), we see that the MC is satisfied for < a < 2. For the ballistic 
motion, Muiz^q n{z) — az, v = 1, a = 2, and R{t — > oo) = (1 + a)~^. This 
system never thcrmalizes to the reservoir temperature, unless it already starts 
at equilibrium. The system acquires an effective temperature different from 
that of the reservoir. This effective temperature is a signature of metastability 
found in glasses, where the FDT does not hold [1,7-9, 12]. 

The first observation of such phenomena was reported by Kauzmann [7]. 
He noticed that when the entropy of a supercooled liquid is extrapolated 
below the glass temperature Tg, it can become smaller than the entropy of 
the cryatalUine solid. To avoid this paradox, he suggested the existence of 
an effective spinodal temperature Tgp in the supercooled liquid phase. Ricci- 
Tersenghi et al. [9] and Cavagna et al. [103] performed single-spin-flip Monte 
Carlo simulations in square lattices with frustration, in which they obtained 
effective temperatures T^f f ^T. Methods for measuring those effective tem- 
perature [5, 10], and many attempts to get a form of FDT for inhomogeneous 
systems have been discussed in the literature [31, 104, 105]. 

It has been shown that a drastic elimination of the fast degrees of free- 
dom in the dynamics of a system may lead to a violation of the fluctuation- 
dissipation theorem [106]. This is due to the fact that equilibration in the 
coarsened description does not necessarily imply full equilibration of the sys- 
tem; therefore a fluctuation-dissipation relation, whose validity is limited to 
equilibrium or local equilibrium states [107-109], may not exist. The theorem 
is valid when there is a great disparity between slow and fast scales in such 
a way that faster scales relax practically immediately. This feature has been 
found in very different situations as in the diffusion of a Brownian particle 
in a shear flow [8,23,110], in the anomalous diffusion problem [1,12,48], in 
systems undergoing activated dynamics [111,112], and in slow relaxation of 
supercooled colloidal systems [113]. This common scenario may suggest that 
the violation of the fluctuation-dissipation theorem could originate from the 
lack of ergodicity inherent to a coarsened description, which is related to the 
tacit reduction of the dimensionality of the system phase space. 
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9 Ballistic Motion 



As discussed earlier in Sec. 4, most of the anomalous diffusion is subdiffusive, 
what can also be observed in most conductors [16]. However, very recently 
in the history of conductivity investigations, supcrdiffusivc and even ballistic 
motion have been produced in laboratories. Indeed, we can find reliable reports 
on ballistic conductivity in carbon nanotubes [54,57], in semiconductors [114], 
and in semiconductor supcrlattices with intentional correlated disorder [55,56]. 

For a simple description of ballistic diffusion, we use Eqs. (36) and (38) 
and obtain 



n{t) = ^ 

TT 



sm{ijj2t) sin(a'it) 



t 



t 



(56) 



The Laplace transform of Eq. (56) gives, as z ^ 0, n{z) ^ z. Consequently, 

V = I and a — 2, which is the ballistic limit. If we set 70 = ttuj2/4:, the initial 
temperature To = 0, in Eq. (55), we get the effective temperature as Tgff as 



A* = 



eff 



1 



2wi 



(57) 



Equation (57) has a control parameter wi, which measures the "hole" in the 
density of states, and how far we are from the result predicted by the FDT. 
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Fig. 3. Normalized mean square velocity as a function of time for the memory given 
by Eq. (56). Here (3 = ^^^2/2 and u>2 = 0.5. Each curve corresponds to a different 
value of a) wi = 0; b) ui\ = 0.25; c) uJi = 0.45. The horizontal lines correspond to 
the final average value As. In agreement with the theoretical prediction. As decreases 
as uji grows. 



Now, we examine the case when A{t) = v{t), the particle's velocity, so 
that we obtain (w^(t) >= (w^)eqA(i). We simulate the GLE for a set of 10000 
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particles starting at rest at the origin, using the memory in Eq. (56) with 
a;2 = 0.5 and different values of uji. The results of these simulations arc 
shown in Fig. 3, where we plot {v'^{t)). We used the normalization {v'^)eq = 1, 
so that A(t). Notice that X{t) does not reach a stationary value; 

rather, it oscillates around a final average value A^. This value of should 
be compared with A* obtained from Eq. (57). 




0.0 ^ ' ' ' ' ^ 

0.0 0.1 0.2 0.3 0.4 0.5 

Fig. 4. A* as a function of the parameter wi. The line corresponds to the theoretical 
prediction given by Eq. (57). Each point corresponds to a value of As obtained from 
simulations like those described in Fig. 3. 



In Fig. 4, we plot A* as a function of ivi as in Eq. (57) with a fixed value 
of a;2 = 0.5. We also plot the average values Ag obtained from simulations for 
different values of oji. Notice that as lui increases, A* decreases as expected. 
The agreement between simulations and Eq. (57) shows that we can predict 
the average value Ag, even when the FDT does not hold. 



10 Shape of things to come 

After the description of some "well-established physics" , we shall take this 
section to discuss some less conventional ideas about irreversibility, diffusion, 
fluctuations, and the approach to equilibrium. 

Blasone et al. [115] have shown that a quantum harmonic oscillator can 
emerge from a couple of classical harmonic oscillators. Every classical oscil- 
lator will obey the laws of classical physics, however, together their behavior 
will follow quantum mechanics. Biro et al. [116], studying the quantization 
of classical fields, demonstrated that a classical system that operates in 5 
dimensions can transmute into a quantum system in 4 dimensions. 

In the same context, the method of Lie symmetry applied to differential 
equations has been often invoked as a mechanism to derive, not only solutions, 
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but also classes of Fokker-Planck equations with non-trivial drift and diffu- 
sion terms [117-125]. The study of Lie-group representations in space-time has 
been recently developed, following similar procedures as those of field theo- 
ries [126]. As a result, the usual Fokker-Planck equation has been derived from 
a U (1) gauge invariant Lagrangian, and the generalization of such a formalism 
for the SU{2) symmetry has provided a new class of Fokker-Planck dynamics, 
which is non-abelian gauge invariant. The drift and the diffusion terms, in this 
situation, are associated with a tensor metric in a Ricmannian manifold. This 
manifold is based on a Galilean metric space-time, introduced (via the light- 
cone of a (4-1-1) Minkowski-like space) to derive the non-relativistic physics in 
a covariant fashion [127 129]. 

The apparent contradiction between the irreversibility of the macroscopic 
phenomena and the time-reversal symmetry of the fundamental laws has 
driven passionate discussion since Boltzmann times. This paradox of irre- 
versibility, discussed in Sec. 7, finds its place both in classrooms and in highly 
specialized conferences. According to Zwanzig [41], there is no paradox. Ac- 
cording to Chaves et al. [130], the paradox still remains as a problem far from 
being solved. Though some people believe that the only mystery related to 
irreversibility is the fact that the universe started in a very special initial 
state, the question is in fact much subtler and deep. The point is that the 
laws of quantum mechanics warrant that time evolution of an isolated system 
is described by a unitary operator that keeps constant the value of the en- 
tropy. The only escape from that fate could be quantum gravity, a theory still 
to be constructed. The point is that the gravitational interaction has infinite 
quantities that have not been renormalized. Thus, we cannot assure that the 
zero-point quantum fluctuations of the gravitational field - which in fact are 
fluctuations of space-time itself! - result only in the renormalization of the 
physically observed quantities. 

Those fluctuations can in fact create a non-unitary contribution to the 
quantum mechanics evolution operator and thus be a fundamental source of 
irreversibility. Chaves et al. [130] suggested that this in fact occurs. On the 
basis of the quantum fluctuations of the metric tensor, they proposed an ex- 
tra term in the Schrodinger equation which makes time evolution operator 
non-unitary. Their calculations demonstrated that the coherence time of the 
microscopic system would be too long to be observed, but macroscopic sys- 
tems would decohere very quickly. Acebal et al. [131] demonstrated that those 
metric fluctuations could also remove the infinities that plague quantum field 
theories. 

These unconventional analyses may prove useful for other derivations. 
They require a deeper understanding of time and space in the field of sta- 
tistical mechanics. 
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11 Spatio-temporal conjecture for disordered system 



Up to now, we have discussed stochastic systems, i.e., systems with noise that 
we shall name temporal disorder. For those, given the NDS, p{uj), it is possible 
to obtain the memory and then, by using Eq. (26), the diffusive exponent a. 

A second class of systems is composed by those which present spatial dis- 
order. They have been thoroughly investigated in the last half century [132], 
nevertheless, some questions concerning localization or diffusion still remain 
open. Let us consider, as an example, the Heisenberg chain [133, 134] 



where S* = 1/2. Here, J; is the exchange integral at the site I. Equivalently, 
we could consider the disordered harmonic chain [135] or even the Anderson 
model [136]. 

Can we predict the properties of those systems in the same way we do for 
the GLE? The answer is partially yes, partially no. The conjecture [84], being 
valid, will help to answer those questions. 

Consider a system which presents fluctuations in its energy density of 
states D[E); let us call pf{E) the fluctuation density, then 



If this is true, then pp can be introduced in Eq. (36) to obtain the diffusive 
exponent a. This is the spatio-temporal conjecture [84]. This conjecture has 
been verified for the quantum disordered Heisenberg chain [137] and is under 
consideration for many similar systems. 

12 Conclusion 

In this review we discussed some old and permanent problems of statistical 
mechanics. The way a system approaches equilibrium is connected with some 
basic questions in physics such as the reversibility paradox, the mixing condi- 
tion, the ergodic hypothesis and the fluctuation-dissipation theorem. We have 
drawn a line between them and, in particular, we have discussed the hierarchy 
MC ^ EH ^ FDT, established by Costa et al. [1]. 

We have approached the problem taking diffusion as a main phenomena in 
physics, since most processes are related with transport of matter, energy, or 
information. In this context, the validity of the FDT is exhibited for ballistic 
motion. Ballistic motion is presented here as the frontier between a stochastic 
process described by a GLE and other processes such as hydrodynamical ones. 

We discussed relaxation processes and the conditions that the correlation 
functions must fulfill. We presented, in Sec. 10, discussions on the frontier 



N 




(58) 



PFiE)=piE). 



(59) 
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of physics with particular consequence to statistical physics. We revive the 
reversibility paradox as an unclosed subject, as well as the MC, EH and FDT. 
We discussed a conjecture that, if valid, will make an important connection 
between stochastic and Hamiltonian descriptions. 

Nonlinear dynamics is a field which deserves much attention; in particular, 
the coalescence of trajectories has been intensively studied in the last few 
years [138, 139]. There, the restriction of the degrees of freedom may confirm 
as well the hierarchy exposed here. 

We have not focused deeper on real complex systems; we chose to follow 
easy-to-understand concepts where limits could be analytically obtained. This 
gave us a good framework for analyzing more complex structures. 

We also tried to show that a given result may be obtained through many 
different formalisms. Feynman once said: "A physicist must know at least five 
different ways to obtain a result" [140]. If we consider the FPE, the FFPE, 
and the GLE as alternative approaches, we are close to fulfilling Feynman's 
requirement. It is nice to know that those approaches agree in the main results; 
however, the full picture has not yet been drawn, particularly for anomalous 
diffusion. 

Although anomalous diffusion remains as a surprising phenomenon, we 
hope that this work will help in the centennial effort to understand diffusion 
and the relation between fluctuation and dissipation. A generalization of the 
FDT to include nonlinearities and ballistic motion is necessary, what will 
require a deeper understanding of systems far from equilibrium. 
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